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Abstract. 

Various sophisticated approximation methods exist for the description of quantum 
many body systems. It was realized early that the theoretical description can simplify 
considerably in one dimensional systems and various exact solutions exist. The focus 
in this introductory paper is on fermionic systems and the emergence of the Luttinger 
liquid concept. 
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1. Introduction 

The theoretical description of systems of interacting particles is notoriously difficult 
and only few exact results are available. Especially in the development of the theory of 
continuous phase transitions it became obvious that the number d of spatial dimensions 
plays an important role. While for most experimentally realizable systems one has 
d — 3, it early turned out that for d = 1, i.e., one-dimensional systems the theoretical 
description simplifies considerably and exact results for various types of spin chains can 
be obtained. By a mapping to fermions this implies exact results also for special lattice 
models of interacting "electrons" in one dimension as discussed in section 2. 
In this introductory paper the focus is on (normal) one- dimensional fermionic quantum 
liquids, which are many body systems in which the indistinguishability of the elementary 
constituents is important. These particles can live on a lattice or the continuous line. On 
low energy scales the metallic state is a non-Fermi liquid characterized by a power law 
decay of space-time correlation functions with interaction dependent exponents. The 
name Luttinger liquid (LL) was termed for this behaviour [I]. 

While in the beginning of the theoretical developments the corresponding models were 
considered a mere playground for theoreticians recent developments have made the 
attempt to experimentally verify Luttinger liquid behaviour a florishing field of research 
as shown in this special issue. 

This paper gives a historical account of the emergence of the Luttinger liquid concept. 
The typical properties were first found in the models proposed by Tomonaga [2] 
and Luttinger [3], which use rather restrictive assumptions about the interaction. 
Tomonaga's important step was to realize and use the fact that the low energy spectrum 
of noninteracting fermions in one dimension is identical to that of a harmonic chain. 
This allows to describe the interacting fermions as a system of coupled oscillators [2]. 
Luttinger's calculation of the momentum distribution in the groundstate marks the 
appearance of power laws for interacting fermions in one dimension [3]. It was realized 
much later that the low energy physics of these models is generic under rather weak 
assumptions [lj H]. This low energy physics can be found also in bosonic many body 
systems [U [5]. Here the focus is on the analytical description of fermionic systems. 
Important computational techniques for one-dimensional quantum many body systems 
like the density matrix renormalization group (DMRG) [6], [7] are not discussed here. 

2. Models 

In this section we present models which played an important role for the theoretical 
understanding of interacting quantum systems in one dimension. We begin with the 
anisotropic spin 1/2 chain with nearest neighbour interaction in an external field h. 
The Hamiltonian reads 

H = £ [J x s*s* +1 + Jys\s y l+l + - hsl] , (1) 
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with the exchange couplings J a . The operators of the spin components on the same 
site i obey the usual angular momentum commutation relations and spin operators on 
different sites commute. For J x = J y = one obtains the Ising chain [8], the simplest 
model of interacting spins. As all operators commute it can be considered a classical spin 
model. Ising's exact solution for the free energy showed that at finite temperatures no 
symmetry breaking to a state with a finite magnetization occurs for vanishing external 
field. It was suspected that this a special property of one dimension [9] , later confirmed 
by Onsager's exact solution for the free energy of the £u>o- dimensional Ising model with 
h = 0, showing a phase transition at finite temperature |10j . 

For the isotropic case J x = J y = J z = J and putting h = one obtains the Hamiltonian 
of the Heisenberg spin 1 /2 chain [11] . For J < the ground state is ferromagnetically 
ordered, while for J > the spins are antiferromagnetically correlated but not ordered 
in the groundstate. In 1931 Bethe [12] presented his famous Ansatz for the exact 
eigenstates and eigenvalues of this model which was later generalized to a larger class 
of l<i-models as discussed in various textbooks [TBI HI IE]- At the end of his paper 
Bethe announced a follow up paper whith a generalization of his Ansatz to higher 
dimensions. The fact that it never appeared, again shows that many body physics in 
one dimension is special. In the following the sophisticated Bethe-Ansatz technique is 
not described, only exact results for Luttinger liquid parameters for lattice models are 
mentioned later. Unfortunately no exact results for the correlation functions discussed 
later are available within the Bethe Ansatz approach, but it should be mentioned that 
significant progress has been made recently to calculate dynamical spin structure factors 
using a non-perturbative ( "form factor" ) approach [T5J [16j [T7] . 

Another special case of the general spin model in Eq. (CD) should be mentioned. For 
Jy — Jz — one obtains the transverse Ising model [18] which now serves as a standard 
model for a system with a quantum phase transition [19]. It was solved exactly [18] 
by the use of a Jordan- Wigner transformation [201 El] which relates the set of spin 1/2 
operators to a set of spinless Fermi operators. For J x = J y the spin model in Eq. (TjQ) 
reads in the fermionic representation 



e o(4 c i - ~ *(4 c i+i + H - c -) + u ( c l c i - \){c\ + ic i+l - -] 



(2) 



where eo = —h,t = —J x /2, U = J z and c\ is the annihilation (creation) operator of 
a fermion at site i . This Hamiltonian describes spinless fermions on a chain with a 
nearest neighbour "Coulomb interaction" . The spin model with J z = ( "Ay-model" ) 
corresponds to noninteracting fermions and can therefore be solved exactly [2TJ . 
The spinless model in Eq. ([2]) is one of the lattice models which played an important role 
in the emergence of the Luttinger liquid concept. The other one is the one-dimensional 
Hubbard model El [23l IE 



H = E [-*( c I<a+v + H - c -) + Uc l^iA,i c i,i\ > ( 3 ) 
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where a is the spin label and U the on-site Coulomb repulsion. 

In his seminal paper Tomonaga [2] treated interacting fermions on the continuous 
line. Without impurities all fermionic models discussed in this paper can be written in 
the form 

H = Yl ek °k c k + 2 Yl v ^k 2 k 3 kAA^ c k S > ( 4 ) 

fc ^1,^2,^3,^4 

where k is a double index k, a for models including spin like the Hubbard model. For 
the lattice models Eq. ([2]) and Eq. ([3]) the momenta ki are in the first Brillouin zone 
and for continuum models they are on the line extending from — oo to oo. The energy 
dispersion e k and the interaction matrix elements v klk2k3k4 are specified in the following 
sections. 

3. The Tomonaga-Luttinger model 

A decisive step towards an understanding of interacting fermions in one dimension 
beyond perturbation theory was Tomonaga's idea [2] to bosonize the Hamiltonian Eq. 
fl4]) for nonrelativistic particles on a line L with periodic boundary conditions (h = 1) 

e k = k 2 /(2m) , v klk2k3k4 = - k 3 )5 kl+k2M+ki ■ (5) 

Tomonaga studied the case when the Fourier transform of the two-body interaction v(k) 
is nonzero only for \k\ < k c = 2irn c /L, where the cut-off k c is much smaller than the 
the Fermi momentum kp = 2irnp/L. This corresponds to a long range interaction in 
real space. Perturbation theory then indicates that the ground state and low energy 
excited states have negligible admixtures of holes deep in the Fermi sea and particles 
with momenta \k\ — kp ^> k c . This is the motivation for Tomonaga's approximation 
to linearize the dispersion e k in the regions around the two Fermi points ±kp, with 
particle-hole pairs present 

k ~ ±kp : e k = €f ± Vp{k =)= kp), (6) 

with vp = kpjm the Fermi velocity. Tomonaga realized that the Fourier components of 
the operator of the density 

fL/2 

Pn = / P {x)e-^ X dx = C l C n> + W (7) 

J-L/2 

where c\,{c n ,) creates (annihilates) a fermion in the state with momentum k n i = 2iin'/L, 
play a central role not only for the interaction but also for the kinetic energy. His 
important idea was to split p n for momentum transfer \k n \ <C kp into two parts, one 
containing operators of "right movers" i.e. involving fermions near the right Fermi point 
kp with velocity vf and "left movers" involving fermions near — kp with velocity — Vp 

Pn=J2 C n' C n'+n + C n' C n'+n = Pn,+ + Pn,- ■ (8) 

n'>0 n'<0 
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In the subspace with no holes deep in the Fermi sea in which all one-particle states \kj) 
with \j\ < M = np — 7n c are occupied, the commutation relations [2] 

[Pm,a, PnA = OLTn8 a p5 m _ n l (9) 

hold for |n|, |m| < M. The dimensionless constant 7 has to be chosen larger for stronger 
interaction. If one defines the operators 

forn > 
for n < 

and the corresponding adjoint operators b\ this leads using pt Q = p_„ iQ , to the bosonic 
commutation relations 

[&«,&J = 0, [6 n ,6] n ]=<5 mn i. (11) 

The kinetic energy of the right movers as well as that of the left movers can be 
expressed as a bilinear form of the 6^ -operators using a remarkable operator identity 
first presented by Kronig in a different context [2U [251 [26] . For the right movers it reads 




27T 

T + = ^v F k n c ] n c n = v F — 

n=l 



^m&J> m + -A/" + (A/"+ + i; 



(12) 



where A/+ = Xl^Li c n c n * s the particle number of the right movers. For its proof the 
commutation relations Eq. ( FlTT) have not to be used. 

As V is bilinear in the p n the same is true for the p nj(Jl . For the linearized fermionic 
dispersion = v F \k\ + const, shown as the dashed line in Fig. 1 and the two-body 
interaction in Eq. (jHJ) the Hamiltonian in Eq. (jlj) can therefore apart from an additional 
term linear in the particle number operators A/± exactly be rewritten as 



n>0 ^ ^ ' 

(b{bl n + 6_ n 6 n ) J + ^ [v N Af 2 + Wjl 7 2 ] = + (13) 



2tt 

where A/" = A/+ + A/1 is the total particle number operator, J7" = A/+ — A/1 the "current 
operator", and the velocities are given by vjy = v F + v(0)/ir and vj = v F . 
If one now assumes that the bosonic commutation relations in Eq. ( 11 ip hold generally 
(see discussion below) the operators H B and Hj^ t j commute and with the Bogoliubov 
transformation cv n = b^ n cosh 9 n — 6_ n sinh 9 n to new boson operators the Hamiltonian 
Hb can be brought into the form 

H B = ^2u(k n )al l a n + const., u(k n ) = v F \k n \y/l + v(k n )/(irv F ) (14) 
and 6 n is determined by 

tanh# n = — v(k n )/(27rv F + v(k n )). (15) 
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Figure 1. Energy dispersion as a function of momentum. The full curve shows the 
usual "nonrelativistic" dispersion and the dashed curve the linearized version. The 
dot-dashed parts are the additional states for fceand = 1.5fci?- The model discussed by 
Luttinger corresponds to fceand — > oo. 



For \k n \ <C k c and a smooth v(k) the boson dispersion is approximately linear u(k n ) « 
^d^nl with v c = y/v^vj the charge velocity. With the approximation to linearize e& 
around the Fermi points it is strictly linear up to k c if v(k) is constant up to the cut-off. 
Besides the charge velocity v c the "stiffness constant" K = \Jvj/vn plays an important 
role. The noninteracting case yields K — 1, attractive interactions v(0) < lead to 
K > 1 and repulsive interaction with v(0) > imply < K < 1. For the generalized 
model, where v(k n ) is replaced by gi{k n ) in the first line on the rhs of Eq. ( [131) (scattering 
events on one of the Fermi points) and by g2{k n ) in the second line (scattering events 
involving both Fermi points which conserve the number of right and left movers) the two 
independent quantities v c and K describe the low energy physics. This generalization to 
non- Galilei-invariant systems turns out to be important for the general Luttinger liquid 
concept discussed in section 4. 

A simple trick to extend the range of validity of Eq. ([111 and therefore for the step 
from Eq. ([131) to Eq. (114j) is to make the Fermi sea deeper for fixed k c by extending the 
left (right) mover branch to positive (negative) fc-values up to ^Band > (— A;B an d < 0) 
as shown in Fig. 1 for the (arbitrary) value fceand = 1.5/cf. The Kronig relation is 
easily extended to this case and leads to an additional term linear in the particle num- 
ber operators and therefore Eq. (fT3|) still holds. Luttinger treated a model with two 
infinite branches of right and left moving fermions with dispersion ±vpk [3]. As he 
made an error related to the fact that his Hamiltonian is not bounded from below, it 
is better to switch from Tomonaga's to Luttinger's model keeping /ceand finite before 
taking the limit /cBand - > [27]. Because of the close relation of both models the term 
" Tomonag a- Luttinger (TL) model" is often used. 



Fortunately Luttinger's error had no influence on his inquiry if a discontinuity of 
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(n k}+ ) a t hp exists in the exact ground state of the interacting model, as expected from 
Fermi liquid theory [281 [29] . After a lengthy calculation using properties of "Toeplitz 
determinants" Luttinger found that the average occupation (71^+) in the ground state 
for k ~ k F in the limit L — > 00 behaves as 

sign(/c F - k), (16) 

where a > depends on the interaction strength (see below). This is the power 
law behaviour mentioned in the introduction. It cannot be obtained by finite order 
perturbation theory in the interaction strength, which produces logarithmic terms in 
I A; — k F \. 

Luttinger's error was corrected by Mattis and Lieb [30] who presented a new algebraic 
method to calculate (n k>+ ). They showed that Eq. (TIB"]) only holds for a < 1. For a > 1 
a term linear in k — k F dominates (n kj +) — 1/2. They also pointed out that one obtains 
(n k)+ ) = 1/2 in the limit interaction cutoff k c — > 00 for a /c-independent interaction. In 
this limit Luttinger's model is equivalent to the massless Thirring model [31] and can 
be written as a (quadratic) local bosonic field theory [HH], not discussed here further. 

Additional information about the system is encoded in its time dependent 
correlation functions (A(t)B), where (. . .) denotes the expectation value in the ground 
state (or in thermal equilibrium) and A(t) = e iHt Ae~ lHt is the operator in the Heisenberg 
picture. As Eq. ( fl4l) implies a n (t) = e~ tu> ^ kn ^ t a n the correlation function (p n (£)/)_.„) can 
easily be calculated using the inverse Bogoliubov transformation. Apart from a prefactor 
its Fourier transform in time is the dynamical structure factor S(q, u>) [32] at q — k n . For 
the Hamiltonian Eq. (114j) S(q, u) is proportional to \q\8(ou—ou(q)). For the nonrelativistic 
dispersion e k = k 2 /(2m) this can only be asymptotically correct for q — > 0. This can 
already be seen from the exact calculation of S(q, u) for the noninteracting case. As 
e kp+q - e kp = qv F + q 2 /(2m) and e kp - e kp ^ q = qv F - q 2 /(2m), for < q < k F the 
dynamical structure factor for fixed q as function of frequency takes the form of a narrow 
box centered at v F q of width q 2 /m and height ~ 1/q. 

In the calculation of one-particle Green functions the Heisenberg operator c n (t) enters. 
It cannot be expressed using the "first step" of bosonization introduced by Tomonaga 
[2]. Even time independent expectation values like the momentum distribution (n kt+ ) = 
(cj, + c k) +) require an additional theoretical concept as used by Luttinger [3] and Mattis 
and Lieb [30] . 

The calculation of (n kt+ ) can be further simplified by bosonizing the field operator. 
This concept was introduced by Schotte and Schotte in the context of x-ray absorption 
from a core hole in the presence of a Fermi sea [33] . To the calculation of correlation 
functions of the Tomonaga-Luttinger-model it was first applied by Luther and Peschel 
[3T| . Later subtleties of this second step of bosonization were addressed [U [351 126] - 
In this step the c n ± are not bosonized directly but the field operators if}±(x), e.g. for 



k — k F 
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the right movers 

1 oo 

^) = -fi E e ^ c «>+ - 

n=— oo 

where the limit /cRand -+ oo apparently is performed first. The commutation relations 
of ip+{x) with the boson operators b m and ¥ m imply the form 



v n 
n=l v 

where the Klein operator + (x) lowers the fermion number by one and commutes with 
all boson operators. Its explicit form has been discussed in detail [H [35j [26] . It is not 
presented here, only the fact that + (x) and 0_(x') anticommute is mentioned. 
Using Eq. (fl8l) . the Bogoliubov transformation and the Baker-Hausdorff formula, 
e A+B = e A e B e~^ A ' B } if the operators A and B commute with [A, B], it is straightforward 
to calculate ground state one-particle Green functions like iG+(x, t) = (^+(0, 0)ip + (x, t)) 
which enter the description of photoemission. Using a n \E ) = one obtains [M] 

ie iia G${x, t) = ^-p- exp i ^ - [ e -^x^ n t) + 2 ginh 2 ^ ( cog ^y^t _ ^ j 

±J I ft 

y n=l 

(19) 

where \x is the chemical potential. Putting t — the momentum distribution (rik t +) 
is obtained by Fourier transformation. In the limit L — > oo the 9 n can be expressed 
by a continuous function 6{k n ) and the equal time Green function G < (x,0) decays like 
x -(i+a) { OI x ^> l/k c with a = 2 sinh 2 [6*(0)], in contrast to x~ x in the noninteracting 
case. Therefore a which in Eq. f|T6|) determines the behaviour of +) near the Fermi 
point kF is called the anomalous dimension. It can also be expressed in terms of the 
stiffness constant K 

a = 2 sinh 2 0(0) — (K — 1) 2 /2K. (20) 

For small interaction a is proportional to v(0) 2 . Only for the special case v(0) = 

the anomalous dimension vanishes and (rik,+) has a discontinuity at kF, the hallmark 

of Fermi liquid theory [28l [29] . In the generic interacting case one has a > and 

Luttinger's power law Eq. (fTB]) holds for a < 1 (see Fig. 2). 

At finite temperatures d{nk,+)/dk diverges like T a_1 at k = kp for a < 1. 

From the Fourier transform of ie t ' lt G < (0, t) one obtains the local spectral density p < (u), 

where u is the energy relative to the chemical potential. For T = this leads to a power 

law suppression [3l] 

" <(w )~ e (-^(S)° (21) 

of the spectral weight near the chemical potential where denotes the unit step function. 
This result holds for arbitrary noninteger values of a > [SB]- The energy range over 
which this asymptotic behaviour can be used depends on the functional form of v(k). 
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v 0.5 




-2-10 1 2 

(k-k F )/k c 

Figure 2. Average occupation (rik : +) for different interactions specified by sinh(9(fc): 
The full line corresponds to sinh 2 6(k) — 0.3e~ 2 l fe l/ fcc i.e. the finite value a — 0.6. The 
expectation from Fermi liquid theory with a finite jump at kp is shown as the dashed 
line for sinh 2 9(k) = 0.6(|fc|/fc c )e _2 ' fc '/ fec i.e. vanishing anomalous dimension. 

For finite temperatures the local spectral weight is nonvanishing also for u > and 
p < (0) ~ T a in the low temperature limit [37] . 

The spectral function p^(k,u) relevant for describing angular resolved photoemission is 
obtained from Eq. (fl9|) by a double Fourier transform. Also the total spectral function 
p + (k,u) = p+(k,u) + p+(k,u) can be obtained from p^(k,u) by using the relation 
P+(kp + k, oj) = p+{k F — k, —oj). 

For a general k- dependence of v(k) the double transform can be performed analytically 
only approximately [36]. The exact calculation is possible numerically, e.g. recursively 
[38| |39| RED"] . An exception is provided by p+(kp,oj). At the Fermi momentum one 
obtains for oj < asymptotically p+(kp,oj) ~ a(— i.e. for 5(0) ^ a power law 
divergence as long as a < 1. There is no sharp quasiparticle peak as in a Fermi liquid. 
This behaviour as well as the power laws in Eq. ( TTBl) and Eq. ( |2TT) are the hallmarks of 
Luttinger liquid behaviour pQ. 

For kp — k c < k < kp and a constant v(q) up to the cutoff fc c the fc-resolved spectral 
function shows a power law singularity at u — v c {k — kp) if a < 1/2 [3H] 

p+(k F + k,u) ~ e(-w - f c |fc|)(-o; + u c fc)f _1 (-a; - v c ^)§. (22) 

This result was first obtained assuming 2 sinh 2 9(q) = ae" r ' <? ' and uj(q) = v c \q\ [31]. A 
comparison of Eqs. (Tl4l) and (IT5|) shows that this is not consistent for all q. As the 
asymptotic analysis in two variables is less developed than for the one-dimensional case 
it is not even known if for k < a power law at oj — v c k exists if the derivatives of v(k) 
are different from zero at k = [36]. Another modification of the power law singularities 
in Eq. ( 1221) for k ^ kp results from the corrections to the linearization of the dispersion 
efc around the Fermi points [HJ H2] as further discussed in section 5. 
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The TL model including spin was introduced and solved exactly using fermionic 
many-body techniques by Dzyaloshinski and Larkin [33] . This included the one-particle 
Green function G±(x,t). The implications for the spectral functions p±(k,u) when spin 
is included were not discussed. 

Using bosonization the step from the spinless to the spinful model is simple. All one has 
to do is to switch from the boson operators bm,a for the two spin components a = j, J, to 
"charge" (c) and "spin" (s) bosons [H H] 

b n<c = (6 n>t + b n ,i) , 6„, s = (6 nt - 6 re> 4) . (23) 

One can write the TL-Hamiltonian H^f^ 1 for spin one-half fermions as [U S] 

H£l 2) = Htl,c + H TL , S , (24) 

where the Htl,ci are of the form of Eq. (1131) but the interaction matrix elements 
have the additional label a. The two terms on the rhs of Eq. (T2"4"|) commute, i.e. 
the charge and spin excitations are completely independent. This is usually called 
spin-charge separation and is another hallmark of LL physics. The "diagonalization" 
of the two separate parts proceeds exactly as before and the low energy excitations 
are "massless bosons" uj nya ~ w a |^n| with the charge velocity v c = [vj c v^ c ) 1 / 2 and the 
spin velocity v s = (f j s f tvJ 1 ^ 2 - The corresponding two stiffness constants are given 
by K c = (vjc/vnJ 1 / 2 and K s = (vjJvn s ) 1 ^ 2 - The low temperature thermodynamic 
properties of the TL-model including spin can be expressed in terms of the four 
velocities vn c , v j c , vn s , vj s or the four quantities v c , K c , v s , K s . For spin rotation invariant 
interactions K s = 1 holds [H H]. 

The one-particle Green functions of the spinful model are given by the square root of 
the product of the charge and spin part which individually are of the form in Eq. (|T9|) 
[1U HU US] • The anomalous dimension is given by a = sinh 2 6* c (0)-|-sinh 2 9 S (0) = a c + a s , 
where a s vanishes in the spin rotation invariant case. Again the spectral function 
p + (k,u) can be calculated analytically in the low energy regime for a constant v(k) up 
to the cutoff k c . It shows two power law singularities jHJ H5] for sufficiently small values 
of the a a . For a s = the "spin singularity" is determined by the exponent (2a — l)/2 
and the "charge singularity" by (a — l)/2. These "peaks" disperse linearly with k — kp- 
For the modification of the singularities due to the k- dependence of v(k) and the correc- 
tions to the linearisation of e& the same arguments hold as in the spinless case. At the 
Fermi momentum one again obtains p + (kp,uj) ~ l^l"" 1 . For the local spectral density 
Eq. ( 12T|) holds also in the spinful model, and the rhs of Eq. ( TT6]) also holds for (n fccr + ). 

Not all interesting results for correlation functions of the TL-model can be listed 
here. The unusual effect of impurites on Luttinger liquids can e.g. be traced back to the 
IQp-^ -1 ) divergence of the static density response function at k = ±2fcp + Q for repul- 
sive interactions [31]. This leads to the breakdown of a perturbational analysis for an 
impurity potential with a weak ±2kp backscattering contribution. The renormalization 
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group analysis by Kane and Fisher [46] showed that the backscattering potential is a 
relevant perturbation for repulsive interactions as expected from earlier work by Mattis 
[17] . The flow to strong coupling implies that the system behaves as if it is split by 
the impurity into two chains with fixed boundaries at the end. Therefore it is necessary 
to mention the behaviour of the one-particle Green function close to a boundary. The 
bosonization for periodic boundary conditions described above has to be modified [48] to 
describe fixed boundary conditions. For spinless fermions and x close to the boundary 
{ip(x, 0)ip^(x, t)) decays like (l/t) 1+olB in the long time limit, where ag = 1/K — 1 is the 
boundary exponent. The local spectral function close to the boundary shows a power 
law p(x,u) ~ M as , where the proportionality factor contains an oscillatory part in the 
position variable. The fact that in the low temperature limit the linear conductance of 
a backscattering impurity vanishes like T 2 ° B [46] can be understood as an end-to-end 
tunneling between the "split chains" . In contrast to the "bulk" anomalous dimension a 
the boundary value is proportional to v(0) for small interactions. 

A challenging problem is to describe the disorder in a Luttinger liquid with a finite 
impurity density [49j [50], [51] . We refer to chapter 9 of Giamarchi's book for an extended 
discussion [4]. 

4. The Luttinger liquid concept 

Tomonaga was well aware of the limitations of his approach for more generic two-body 
interactions ("In the case of force of too short range this method fails" [2]). It was only 
realized later that the TL model is the fixed point Hamiltonian for a rather general 
class of models [521 [TJ [531 El]- This emergence of the general Luttinger liquid concept 
is discussed in this section. 

In the opposite limit k c ^> hp and a ^-independent interaction Tomanaga's continuum 
model corresponds to a short range interaction in real space. Then the low energy 
scattering processes with momentum transfer ±2/cp have to be included. They are 
usually modeled by the additional u gi' -interaction term 



Introducing a band cutoff Solyom [52] made a renormalization group (RG) study of 

this interaction at the one loop level. If the variable s runs from zero to infinity in 
the process of integrating out degrees of freedom he obtained for spin-independent 
interactions gq = g iX = & (i = 1, 2) 




(25) 




dgi(s) 




dg 2 {s) 




(26) 



ds 



ds 



with the solution g\{s) = g\/[l + sgi/ (nvp)}, where g x is the starting value. The 
(^-interaction is not renormalized. For gi > the interactions flow to the fixed line 
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gl = 0,g% = g2 — gi/2 and the fixed point Hamiltonian is a TL-model [22]. This shows 
the generic importance of the TL-model for repulsive interactions. 

For #1 < the flow is to strong coupling. In order to understand the strong coupling 
regime Luther and Emery bosonized the additional interaction in Eq. ( 125]) and 
showed that "spin-charge separation" also holds for this model [55]. The charge part 
stays trivial with massless charge bosons as the elementary interactions. They showed 
that for a particular value of g\\\ the exact solution for the spin part of the Hamiltonian 
is possible using refermionization. The spectrum for the spin excitations is gapped. It 
is generally believed that these properties of Luther- Emery phases are not restricted to 
the solvable parameter values. 

Strong coupling phenomena which lead to deviations from LL behaviour can occur 
in the lattice models like the ones discussed in section 2, when for commensurate fillings 
Umklapp processes become important [I]. Here important results for the models in Eqs. 
(EJ) and ([3]) are presented. 

For the spinless fermions Eq. ([2]) the parameters in Eq. ((4]) for a chain of N sites 
with periodic boundary conditions and fc-values in the first Brilloin zone are given by 



The m = term on the rhs of Eq. (j27|) represents the direct scattering terms and 
the m — ±1 terms the Umklapp processes. It is a low energy process in the half 
filled band case hp = it/2 discussed here. Renormalization group analysis around the 
noninteracting fixed point shows that the Umklapp terms are strongly irrelevant which 
implies that the system is a LL for small U > [56]. For U 3> t > charge density 
wave (CDW) order develops in which every other site is occupied in order to avoid the 
Coulomb penalty. The mapping to the spin model suggests that the transition occurs 
at U c = It as for U c > 2t the Ising term dominates. The exact Bethe ansatz solution 
confirms this and shows that the model at half filling is a LL for \U\ < 2t. The Luttinger 
liquid parameters can be obtained using a ground state property and the lowest charge 
excitation [57]. For repulsive interactions U > the value of K = 7r/[2arccos (— U/2t)] 
decreases monotonously from the noninteracting value K = 1 to K = 1/2 for U = 2t, 
which corresponds to an anomalous dimension a = 1/4. In order to reach smaller values 
than 1/2 for K the interaction has to have a longer range in real space [I]. The limit 
K — > is reached by the bare Coulomb interaction as v(k) ~ log (l/|/c|) for k — > and 
the system is not a LL. The 4/c^-harmonic of the density- density correlation function 
shows a very slow decay almost like in a Wigner crystal [58] . 

The limit in which the lattice constant and the density go to zero corresponds to the 
continuum limit. The interaction goes over to a contact interaction. Because of the 
Pauli principle its effect vanishes and K — > 1. This limit is very different for the Hub- 
bard model Eq. fl3]) as the onsite interaction is between electrons with different spins. 



6k = —2t cos k , 




m=0,±l 



(27) 
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The energy dispersion e& for the Hubbard model Eq. ([3]) is the same as in Eq. 
( 1271) and the interaction matrix elements vu^m^M^M^a have the same m-sum. The 
fc-independent prefactor is proportional to S^^S^^d^^^ . To show the difference to 
the spinless model the focus is again on the half filled band case which is metallic for 
{7 = 0. The limit U ^> t is easy to understand. Each site is singly occupied in order 
to avoid the Coulomb penalty. In this limit the model can be mapped to the spin 
model Eq. (CQ) with J x — J y — J z — 4t 2 /U and h = 0, i.e. the spin-1/2 Heisenberg 
antiferromagnet which has gapless excitations [131 IH E] • I n contrast there is a large 
gap A c « U for excitations in the charge sector. As the model can be solved exactly by 
a generalized Bethe ansatz approach this Mott-Hubbard gap can be obtained exactly by 
solving Lieb and Wu's integral equations [23j[14]. It turns out to be finite for all U > 0. 
It is exponentially small A c m (8t/ir)^U/texp (— 2iit/U) for < U <Ct. This shows 
that the Umklapp term is no longer irrelevant at the noninteracting fixed point. As the 
Pauli principle does not influence electrons of opposite spin the RG analysis shows that 
the Umklapp terms are marginally relevant at the noninteracting fixed point [591 111 EH] ■ 
When the band is not half filled Umklapp is not a low energy process and the Hubbard 
model is a Luttinger liquid with K s = 1 for all U > 0. The LL parameters K c and 
v a (a = c, s) can be obtained by numerically solving Lieb and Wu's integral equations. 
The results show that K c — > 1/2 for n — > as well as n — > 1 (half filling) for all U > 
[601159] . 

For results for various correlation functions we refer to the textbook on the one- 
dimensional Hubbard model |14j . 

The discussion of the two lattice models Eqs. ([2]) and ((3|) shows the general im- 
portance of the Luttinger liquid concept for one-dimensional fermions with repulsive 
interaction. Only for half filling qualitative deviations occur. 

As mentioned in the introduction the focus of this paper is on fermionic systems. 
Bosons in one dimension with repulsive interaction also behave as Luttinger liquids. 
For them there is an essential difference between the noninteracting and the interacting 
system [HEl]. For the nonrelativistic dispersion = k 2 /(2m) the excitation spectrum 
of the bosonic many body system is linear in \k\, i.e. has the typical LL form only if a 
finite interaction is present. In addition to the sound velocity the low energy physics is 
described by the stiffness constant K which again determines the large distance and 
long time behaviour of correlation functions. In contrast to fermionic systems the 
noninteracting bosons correspond to the limit K — > oo. 

Very versatile systems to experimentally test the Luttinger liquid behaviour for bosons 
in one dimension are ultracold gases in strong optical lattices [5]. 

Examples for the experimental realization of Luttinger liquid behaviour in quasi- 
one-dimensional electronic sytems are presented in a separate introductory paper 
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[52] . The comparison of theory and experiment faces the problem that strictly one- 
dimensional systems are a theoretical idealization. Apart from this even the coupling 
to an experimental probe presents a nontrivial disturbance of a Luttinger liquid. The 
coupling between the chains in a very anisotropic 3c? compound generally, at low enough 
temperatures, leads to true long range order. The order develops in the phase for which 
the algebraic decay of the corresponding correlation function of the single chain LL is 
slowest [59] . This can lead e.g. to charge density wave (CDW), spin density wave (SDW) 
order or superconductivity. Unfortunately the weak coupling between several LLs or the 
coupling of a LL to a substrate is theoretically not very well understood [H [26]. The 
discussion could easily fill a paper itself. 

Carbon nanotubes are one-dimensional metallic systems when the semimetallic car- 
bon sheet is properly wrapped. In contrast to the systems discussed above two bands 
cross the Fermi level like in a two-leg-ladder [I]. Within an approximate treatment of 
the interaction terms it is possible to describe the system in close analogy to the "reg- 
ular" Luttinger liquids discussed so far [631 EI]- 

One-dimensional metals having elementary excitations which propagate along the 
boundary of a two-dimensional system in one direction only are another type of Luttinger 
liquids. Wen introduced the concept of chiral Luttinger liquids to describe the edge 
excitations in the fractional quantum Hall states [65] [66]. 

5. Extensions of the Luttinger liquid concept and outlook 

The Luttinger liquid concept has recently been extended further in various ways. An 
important step was to examine effects beyond Tomonaga's linearization Eq. ([6]) by in- 
cluding terms of order (k — kp) 2 (or order (k — kp) 3 for the half filled lattice models) 
[4T1 14*2]. This makes relaxation processes possible which do not exist in "linear" Lut- 
tinger liquids. The interaction modifies e.g. the "narrow box" centered around qvp in 
S(q,oj) discussed in sec. 3 and the fc-resolved spectral functions p(k,cu) for k ^ kp. The 
new "nonlinear Luttinger liquid" phenomenology makes contact to methods developed 
for describing the x-ray edge singularity of core hole spectra in the presence of a Fermi 
sea [6T[ 168] . The many-body dynamics is described using effective models for mobile 
quantum impurities in a linear Luttinger liquid [691 HD 112] • 

The temperature dependence of the spectral functions was discussed only briefly so 
far. In the low temperature regime ksT <C v c k c the power law behaviour is smoothed 
out but the anomalous dimension can be recovered using the low energy scaling relation 
P < (oj, ksT) = T a F(u /ksT) [37], [TO], [Tl] . A new scenario can result in the limit that one 
of several intrinsic energy scales goes to zero. An example is the half filled Hubbard 
model discussed in the previous section. In the limit U — > oo the exchange coupling 
J = At 2 /U goes to zero and with it the spin velocity v s . In the temperature range 
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J <C ksT <C v c /a called spin-incoherent Luttinger liquid regime the one-particle 
spectral functions show qualitatively different behaviour from that discussed in section 
3 [?2l Details can be found in a review by Fiete [73] . 

Helical Luttinger liquids are realized in helical conductors e.g. on the edges of topological 
insulators. Unlike the chiral Luttinger liquid a helical LL does not break time reversal 
symmetry in these systems with strong spin-orbit interaction. Helical LL exhibit spin- 
filtered transport where right movers carry spin up and left movers spin down [75] [76] . 
Yet another type of Luttinger liquid behaviour can be realized in the presence of nuclear 
moments. A new ordered phase can result by the coupling of these moments to the con- 
duction electrons via the RKKY interaction [77]. The resulting low energy physics was 
dubbed spiral Luttinger liquid by the authors. A detailed comparison of the spectral 
properties of regular, helical and spiral Luttinger liquids was presented recently [711 [79] . 

In recent years a very active field of theoretical research is to generalize the descrip- 
tion of one-dimensional quantum systems to conditions far from thermal equilibrium. 
A typical example is a finite (interacting) quantum wire which in the initial state is 
attached from the left and right to noninteracting leads with differing chemical poten- 
tials fiL(R) and temperatures T L ^ [801 EE] • The theoretical description is usually done 
using the Keldysh technique [821 [83] . For an even more general class of nonequilibrium 
states where the initial states of the leads are not of the grand canonical form a new 
bosonization technique has been developed [S3] which uses concepts familiar from "full 
counting statistics" [H5] as well as the x-ray edge problem [671 EE] . 

Obviously this short introduction cannot cover all the important contributions 
to the theory of one-dimensional quantum many body systems. As we started the 
discussion with the spin 1/2 chain we end by mentioning that important insights into 
the physics of isotropic antiferromagnetic chains of arbitrary spin were obtained using 
methods of conformal field theory [861 EZl [88] . 
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